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Abstract. We propose a random matrix model for families of elliptic curve 
l^ . L-functions of finite conductor. A repulsion of the critical zeros of these L- 

Cn ' functions away from the center of the critical strip was observed numerically 

C_) . by S. J. Miller in 2006 fSO"; such behaviour deviates qualitatively from the con- 

j^ ' jectural limiting distribution of the zeros (for large conductors this distribution 

)-^ . is expected to approach the one-level density of eigenvalues of orthogonal ma- 

r~^ ' trices after appropriate rescaling). Our purpose here is to provide a random 

matrix model for Miller's surprising discovery. We consider the family of even 
quadratic twists of a given elliptic curve. The main ingredient in our model is a 
calculation of the eigenvalue distribution of random orthogonal matrices whose 
characteristic polynomials are larger than some given value at the symmetry 
r~| . point in the spectra. We call this sub-ensemble of S0(2A^) the excised orthogo- 

nal ensemble. The sieving-ofF of matrices with small values of the characteristic 
polynomial is akin to the discretization of the central values of L-functions im- 
plied by the formulae of Waldspurger and Kohnen-Zagier. The cut-off scale 
appropriate to modeling elliptic curve L-functions is exponentially small rela- 
Cp ' tive to the matrix size N . The one- level density of the excised ensemble can 

be expressed in terms of that of the well-known Jacobi ensemble, enabling the 
former to be explicitly calculated. It exhibits an exponentially small (on the 
t:;J^ . scale of the mean spacing) hard gap determined by the cut-off value, followed by 

■^ ' soft repulsion on a much larger scale. Neither of these features is present in the 

r — I , one- level density of SO{2N). When iV — >■ oo we recover the limiting orthogonal 

f^ ' behaviour. Our results agree qualitatively with Miller's discrepancy. Choosing 

the cut-off appropriately gives a model in good quantitative agreement with the 
number-theoretical data. 
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1. Introduction 

Understanding the ranks of elliptic curves is a central problem in number the- 
ory. The celebrated Birch and Swinnerton-Dyer conjecture provides an analytic 
approach to studying them via the order of vanishing of the associated L-functions 
at the centre of the critical strip: it states that the order of vanishing of the el- 
liptic curve L-function at the central point equals the rank of the Mordell-Weil 
group. On the other hand, following the Katz-Sarnak philosophy, zero statistics 
of families of L-functions are conjectured to have an underlying symmetry type 
associated with certain random matrix ensembles [IHlllI]: in an appropriate as- 
ymptotic limit, the statistics of the zeros of families of L-functions are modeled 
by (a subgroup of) one of the classical compact groups U(A^), O(A^) or USp(2A^). 

If we could model the order of vanishing for families of elliptic curve L-functions 
in terms of random matrices, then we would also gain information about the distri- 
bution of ranks within a family of elliptic curves via the Birch and Swinnerton-Dyer 
conjecture. This is the wider goal that motivates the current work. 

There has been much success modelling families of L-functions with matrices 
from the classical compact groups [71[8l[l0Hl9l[2ll[2?,'25i-'35i,'39i-'33i,'37l[5ll5l[5^ 
including families of elliptic curve L-functions [13, 14,37,38,48,49,6^ . This work 
strongly supports the Katz-Sarnak philosophy that, in the correct asymptotic 
limit, matrices from the classical compact groups accurately model the zero sta- 
tistics of L-functions. Specifically, for families of elliptic curve L-functions it has 
been shown for restricted test-functions that the one- and two-level densities do in- 
deed show the expected symmetry of orthogonal typqj, that is, of random matrices 
from O(A^), in the limit of large conductor |18 | H9 | [62]. 

For the Riemann zeta function, and more generally families of L-functions, 
good qualitative agreement is also found, as the asymptotic limit is approached, 
by modeling the L-functions with finite-sized matrices where the matrix size is 
chosen so that the mean density of matrix eigenvalues matches the mean density 
of zeros of the L-functions under consideration (originally proposed by Keating and 
Snaith [12|I^ ). This finite matrix size model has been refined still further with the 
incorporation of number-theoretic information, leading to extremely convincing 
models for, for example, correlation functions [HE], spacing distributions [5], and 
moments and ratios of values of L-functions [m[T2] . 

However, in the elliptic curve case S. J. Miller, in numerical data gathered 
in [50j , discovered a significant discrepancy from the scaling limit of the expected 
model of orthogonal matrices. That is, for elliptic curve L-functions of finite 
conductor there is a large (repulsive) deviation of the zero statistics from the 
expected orthogonal symmetry that is not explained by either using finite-size 
matrices or by the inclusion of arithmetical terms as in [l]46|ITl]IT2] . Miller was the 
first to discover the soft repulsion of zeros of elliptic curve L-functions from the 
central point, and his data is reproduced in figure |2j However, the larger data sets 



If the family has geometric rank r, then by the Birch and Swinnerton-Dyer conjecture and 
Silverman's specialization theorem all but finitely many of the elliptic curve L-functions have 
r zeros at the central point; the correct scaling limit is the limit of block diagonal orthogonal 
matrices with an r x r identity matrix as the upper left block. 
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of the present paper (see figure Hl|, along with the random matrix model, indicate 
a hard gap containing no zeroso The specific goal of this paper is to present a 
model that accurately describes the zero statistics of elliptic curve L-functions of 
finite conductor. 

The main ingredient in the new model is to mimic central values of elliptic curve 
L-functions by the value of the characteristic polynomial at 1 for matrices from 
S0(2iV). The key fact about these central values, not previously incorporated 
into random matrix models for the zero statistics, is that they are either zero 
(for curves with rank greater than zero) or they are greater or equal to a finite 
minimal size fixed by arithmetical considerations. The characteristic polynomial 
of a matrix A G S0(2A^) is given by 

N 

AA(e^^ N) := det(/ - e''A-^) = JJ(1 - e*(^-^^))(l - e"^'+'''>), (1.1) 

fc=i 

with e^*^^ , • • • , e^*^^ the eigenvalues of A. Motivated by the arithmetical size con- 
straint on the central values of the L-functions we seek to understand, we here 
consider the set T;^ of matrices A G S0(2A^) whose characteristic polynomial 
Aa{1,N) is larger than some cut-off value exp{X). We call this the excised or- 
thogonal ensemble. In figure [T] we plot the distribution of the first eigenvalue of 
5*0(24), which shows no repulsion, as well as the distribution of the first eigenvalue 
of an excised ensemble of 24 x 24 orthogonal matrices. Here the size of the excision 
is exp(A:') ^ 0.005 and we see the hard gap at the origin where no eigenvalues lie, 
as well as the soft repulsion that interpolates between zero and the bulk of the 
distribution. This is seen more clearly in the right hand figure where a region near 
the origin is enlarged. 

Recall that the one- level density Ri for a (circular) ensemble G{N) of matri- 
ces whose eigenvalues are parametrized by an unordered A^-tuple of eigenphases 
{6'„}f is given by 

/?f(^)(^) = iv /... f p{e,e2,...,eM)de2...deM (1.2) 

where P{6,62, . . . ,9n) is the joint probability density function of eigenphasecl- 

T' 

Helpfully, the probability density function of the one- level density Ri'^ can be 
expressed in terms of the well-known Jacobi ensemble Jn (see [23] for properties 
of this ensemble): 

Theorem 1.1. The one-level density R^^ for the set Tx of matrices A G S0{2N) 
with characteristic polynomial A^(e*^, A^) satisfying log |A^(1, A^)| > X is given by 

RT.^e^) = ^ r^^ ^;v. exp(-^^) ^J^(g^. ^ _ ^/2, -i/2)dr 
2ni .L_,^ r 



^Subsequent to our posting this manuscript on the arXiv, Marshall "^S] has developed a 
rigorous theory for this hard gap. 

"^Note that in the case of our interest, namely the (full, or excised) ensemble of orthogonal 
matrices of size 2iV, there are N pairs of eigenvalues parametrized by N eigenphases. 



DUENEZ, HUYNH, KEATING, MILLER, AND SNAITH 






Figure 1. The distribution (blue, solid curve) of the first eigen- 
value of 5*0 (24), showing no repulsion, and the distribution (red 
dots) of the first eigenvalue of an excised ensemble of 24 x 24 orthog- 
onal matrices with exp(A') ^ 0.005. The right hand figure shows an 
enlargement near the origin. The 5*0(24) curve is computed using 
the numerical differential equation solver for Painleve VI which is de- 
veloped in |20]. For the excised ensemble, 3 000 000 random 5*0(24) 
matrices were generated and those not satisfying Aa(1, A^) > exp{X) 
were discarded. 

where Cx is a normalization constant defined in Ii6.12\) and 

/?^(^i;r-l/2,-l/2) = iV /""■•■ [^ flw^'-^^^^-^^^^cosej) 

Jo Jo j^i 

X JJ(cos^j-cos^fe)2d^2---rf^iv (1.3) 

is the one-level density for the Jacobi ensemble Jn with weight function 
u;("'^)(cose) = (1 - cos^)"+^/2(l + cose)'^+^/2^ a = r -1/2 and P = -1/2. 

Applying the method of orthogonal polynomials we evaluate /2i^(6'; a, /3) to give 
Theorem 1.2. With R^^ as above, we have 



«-'<^'-i/ 



"+*°°exp(-rA'), 



^N^+2Nt-N 



^^ ^^ r(2 + j)r(i/2 + j)r(r + 1/2 + j) ^^ 

X (l-cos^)"^ = Miv ^1 lyv ^^; P(N,r,e)dr 

^ ' 2N + r-lT{N + r-l/2)T{N -1/2) ^ ' ' ^ 

with normalization constant Cx defined in Ii6.12\) and P{N,r,6) defined in terms 
of Jacobi polynomials in Ii6.33\) . 

An immediate consequence using residue calculus is the following: 
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Theorem 1.3. With R^^ as above, we have 



oT...^ fo fordie,X)<0 

1 y) \Rl''^''^\e) + C;,ZT=obkiO)eMik + l/2)X) fordie,X)>0, 

(1.5) 
where d{6, X) := {2N — 1) log 2 + log(l — cos^) — X and bk are coefficients aris- 
ing from the residues. The normalization constant Cx is defined in \6.1^) and 
Ri (6) is the one-level density for SO {2 N), defined by U.^) . 

Thus in the hmit X ^ -oo, 6 fixed, Ri^'iO) -^ Rl^^'^^\e). 

In order to compare our random matrix resuhs with the elliptic curve data we 
shall need to take —X proportional to A^, i.e. the cut-off is exponentially small 
in N. From Theorem II. 3[ and illustrated in figure [13], one sees the existence of a 
hard gap where R^^ (9) vanishes, namely the interval {^ > | d(9, X) < 0}. When 
~X oc N the hard gap is therefore exponentially small in N . Beyond this range 
the formula exhibits soft repulsion of the eigenvalues from 1 on a much larger 
scale; the repulsion extends far beyond the hard gap. It is interesting that such a 
'pierced' subset of S0(2A^), where the cut-off is exponentially small on the scale 
of the mean eigenvalue spacing, has such a pronounced effect on the eigenvalue 
statistics over a much larger distance. 

This agrees qualitatively with the discrepancy Miller observed. We go on to 
discuss two ways to determine the relevant parameters which lead not just to 
qualitative but also quantitative agreement. 

In the following sections we give some background information on elliptic curves 
(section [2. II) . detail the anomalous zero statistics observed by Miller in [50] (section 
12.21) and set out the new model (section [3]). The model involves selecting from 
S0(2A^) those matrices whose characteristic polynomial at 1 is larger than some 
given cut-off value (section l5.ip . We apply the model with two matrix sizes: one 
the standard value related to the mean density of zeros and one an effective matrix 
size determined from lower-order terms of the one- level density (section H]). We 
present numerical evidence for our model using the distribution of the first zero of 
a family of elliptic curve L-functions (sections 15. 2 1 15. 3 p and also using the one- level 
density statistic for the excised ensemble of matrices (section [6]). 

2. Elliptic curve L-functions 
2.1. Background. An elliptic curve E can be written in Weierstrass form as 

y'^ + Cixy + c^y = x^ + C2X^ + C4X + Cg, Q G Z. (2.1) 

The L-function Le{s) associated with E is given by the Dirichlet series 

i-W = E^. (2-2) 

n=l 

where the coefficients {X{n) = a{n)/y/n, with a{p) = p + 1 — ^E(Fp), ^E(Fp) 
being the number of points on E counted over Fp), have been normalized so that 
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the functional equation relates s to 1 — s: 

M.)^.(E)(J|)^'"Mi^L.(l-.). (2.3) 

Here M is the conductor of the elliptic curve E\ for convenience we will consider 
only prime M. Also, uj{E) is +1 or —1 resulting, respectively, in an even or 
odd functional equation for Le- There is a Generalized Riemann Hypothesis for 
L-functions of elliptic curves, stating that the non-trivial zeros of Le{s) lie on 
the critical line (with real part equal to 1/2). The even (odd) symmetry of the 
functional equation implies an even (odd) symmetry of the zeros around the central 
point s = 1/2, where the critical line crosses the real axis. Often when we consider 
zeros statistics of a family of L-functions we are particularly interested in the zeros 
close to the central point. 

The family of elliptic curve L-functions for which numerical evidence is presented 
in this paper is that of quadratic twists of a fixed curve E. Let LE{s,Xd) denote 
the L-function obtained by twisting Le{s) by a quadratic character Xd- Here d is 
a fundamental discriminant, i.e., d G Z such that p"^ \ d for all odd primes p and 
d = 1 mod 4 or (i = 8, 12 mod 16, and Xd is the Kronecker symbol (an extension 
of the Legendre symbol, taking the values 1, 0, or —1). The twisted L-function, 
which is itself the L-function associated with another elliptic curve E^, is given by 

Lds. X.) = E ^^%^ - n (l - ^%^ + '^■"'rtx^(ri' ) " , (2.4) 

n=l ^ p V P P / 

where V'm is the principal Dirichlet character of modulus M: 

{liip\M 
^m{p) = {r, ,. . (2-5) 

I U otherwise. 

The functional equation of this L-function is, for (c?, M) = 1, 

L.(s, X.) ^ U-MME) (^)"' ^ffiyiji^d - '• X.y (2.6) 

The sign of this functional equation is Xd{—M)uj{E) and it is more instructive 
to restrict to the fundamental discriminants which, for a fixed E, give an even, or 
alternatively an odd, functional equation for Le{s, Xd)- In particular, the families 
we will consider in this paper will be denoted J-'^(X): those quadratic twists 
of the curve E that have an even functional equation and < rf < X (or else 
—X < d < 0). In the appropriate asymptotic limit, zero statistics of families of 
L-functions of elliptic curves with even functional equation are believed to follow 
the distribution laws of eigenvalues of the even orthogonal group S0(2X) and 
those with odd functional equation are expected to show S0(2X -|- 1) statistics^ 
The asymptotic parameter of the family J^^(X) is X. The conductor of Le{s, Xd) 
is Md^, if M is the conductor of the original curve, and so d is the parameter 
that orders the curves by conductor. It is expected that as X — )■ cxd the zero 
statistics of J-'^(X) tend to the large-X limiting statistics of eigenvalues from 



If the family has rank then we must modify the random matrix ensemble as in Footnote [TJ 
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S0(2A^) in accordance with the Katz-Sarnak philosophy. We propose a model for 
the behaviour of zero statistics for finite X. 

Key to this model are formula of Waldspurger jBU], Kohnen-Zagier [H] and 
Baruch-Mao [3] which show that modular L-functions can only attain discrete 
values at the centre of the critical strip. In particular we have for the twists of an 
elliptic curve L-function 

ceMY 



LE^-l'^^Xd) = K.E- 



rfV2 



(2.7) 



where CE(|(i|) are integers and the Fourier coefficients of a weight-3/2 modular 
form and ke is a constant depending on the curve E. 

2.2. Unexpected repulsion. In 2006 S. J. Miller [5U] investigated the statistics 
of the zeros of various families of elliptic curve L-functions. Figure [2] gives an 
example of what he discovered. It shows a histogram of the first zero above 
the central point for rank zero elliptic curve L-functions generated by randomly 
selecting the coefficients ci up to cq (as defined in (12.11) ) for curves with conductors 
in the ranges indicated in the caption. The zeros are scaled by the mean density 
of low zeros of the L-functions and the plots are normalized so that they represent 
the probability density function for the first zero of L-functions from this family. 
Miller observes that there is clear repulsion of the first zero from the central point; 
that is, the plots drop to zero at the origin, indicating a very low probability of 
finding an L-function with a low first zero. 
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Figure 2. First normalized zero above the central point: Left: 750 
rank curves from y'^ + aixy + a^y = x^+a2X^+a4^x+aQ, log(cond) G 
[3.2,12.6], median = 1.00, mean = 1.04, standard deviation about 
the mean = .32. Right: 750 rank curves from y'^+aixy+a^y = x^+ 
a2X^ + a^x + Uq, log(cond) G [12.6, 14.9], median = .85, mean = .88, 
standard deviation about the mean = .27 



What is surprising about these plots is that the standard way to model such 
L-functions would be with matrices from S0(2A^), with N chosen to be equal 
to half the logarithm of the conductor of the curve. This choice of N has the 
effect of equating the density of eigenvalues to the density of zeros near the crit- 
ical point and there has been much work showing, for the Riemann zeta func- 
tion [T5l[IM9l[22l|271[29l|3ll[33lll2lll7] and for families of L-functions not neces- 
sarily associated with elliptic curves [71El[inHIl[ISl[IIlElll31ISB] , that this can be 
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effective at modelling number-theoretic data even far from the asymptotic limit 
specified by Katz and Sarnak. However, as figure [3] illustrates, even for small-size 
S0(2A^) matrices, the distribution of the eigenvalue closest to 1 on the unit circle 
(the random matrix equivalent to the distribution of the first zero of L-functions) 
shows no repulsion at the origin of the distribution — a fact of course well-known 
in random matrix theory. 
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Figure 3. Probability density of normalized eigenvalue closest to 
1 for SO (8) (sohd), SO (6) (dashed) and SO (4) (dot-dashed). 



This, then, is the mystery. No one doubts that in the large-conductor limit the 
distribution of the first zero of elliptic curve L-functions will tend to the large- A^ 
limit of the distribution of the first eigenvalues of S0(2A^) matrices: Miller observes 
(see figure [2]) that the repulsion decreases with increasing conductor. However, for 
finite conductor we do not get qualitative agreement with the statistics of finite- 
sized matrices, as we do in many other families of L-functions. Indeed one sees 
a qualitatively unexpected phenomenon, namely, repulsion of zeros away from 
the central point. In section E] we outline an idea that successfully models the 
anomalous statistical behaviour discovered by Miller. 

In this paper we concentrate on developing a model for the unexpected zero 
statistics of rank zero curves seen in figure |2l However, it should be noted that 
Miller ^50j also investigated zeros of higher rank curves and the behaviour of zeros 
further from the central point, which will guide future attempts to model statistics 
of L-functions associated to elliptic curves of higher rank. Briefly, the numerical 
findings were the following: 
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(1) The repulsion of the low-lying zeros increased with increasing rank, and 
was present even for rank curves. 

(2) As the conductors increased, the repulsion decreased. 

(3) Statistical tests failed to reject the hypothesis that, on average, the first 
three zeros were all repelled equally (i. e., shifted by the same amount). 

3. The model 

We propose to model the zero statistics of rank-zero elliptic curve L-functions 
with the subset T;^ of S0{2N) defined in the introduction. The justification 
for employing an excised random matrix model such as this is that formula fl2.7p 
indicates that the L- functions themselves have a discretization at the central point. 
Thus the statement 

L^(l/2,x.)<^ (3.1) 

implies that 

LE{l/2,Xd) = 0. (3.2) 

Hence, each Le{s, Xd) £ J^e associated with a rank curve Ed satisfies 

L^(l/2,x.)>|^. (3.3) 

Thus, since previous work (TUIIISIEIIIS] shows that values of L-functions at the 
central point can be characterized using characteristic polynomials evaluated at 
the point 1, to model rank curves we discard from our orthogonal ensemble all 
matrices not satisfying (13. 4p and we will propose a value of X based on (13. 3p . For 
a start, if we are working with discriminants of size around d then equating the 
density of eigenvalues, N/n, with that of zeros near the central point, ^ log(^^^^^), 
gives us an equivalent value of A^: A^std ~ log d. Since the L- values are discretized 
on a scale of l/y/d, we "excise" (i. e., discard), characteristic polynomials whose 
value at 1 is of the scale exp(— A^std/2). 

Thus, for certain values of A^ we plan to model elliptic curve zero statistics using 
matrices from S0(2A^) satisfying 

|Aa(1, A^)| > expA' = c X exp{-N,tj2), (3.4) 

and we will propose a value for c in section 15.11 We present the model for the 
zero statistics using two different matrix sizes. One obvious choice is A^ = A'std 
and we present data for this case in section 15.21 The other choice is an "effective" 
matrix size, Aefr, that incorporates arithmetic information and is based on an 
idea of Bogomolny, Bohigas, Leboeuf and Monastra [5]. We illustrate how A'eff is 
calculated in sectionlHand present the numerical results of the model in section [3751 
In this paper we compare this model with numerical data from a family of twists 
of one particular elliptic curve. En. 

4. lower-order terms and the EFFECTIVE MATRIX SIZE 

By the Katz-Sarnak philosophy the statistical properties of zeros of L-functions 
should asymptotically behave like the scaling limits of eigenvalues of random ma- 
trices drawn from one of the classical compact groups. However, for finite values of 
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the asymptotic parameter we observe deviations (of number-theoretic origin) from 
the hmiting result. We note that for many famihes of L-functions these deviations, 
at least over length scales on the order of the mean spacing of zeros, are slight 
perturbations of the limiting results, in contrast to the elliptic curve case where 
we see distinct, qualitative disagreement (the repulsion at the central point) |3 

In this section we adapt the method of Bogomolny, Bohigas, Leboeuf and 
Monastra [5], previously used to investigate statistics of zeros of the Riemann 
zeta function at a finite height on the critical line, and apply it to the family 
J-^{X) of elliptic curve L-functions. 

Specifically, Bogomolny et al. use a conjectured formula for the two-point cor- 
relation function of the Riemann zeros at finite height on the critical line [6J and 
compare it to the finite- A^ form of the two point correlation function of unitary 
matrices from U(A^). Under the standard procedure of equating mean densities of 
eigenvalues with the mean density of zeros, yielding in this case 

iV = log(^). (4.1) 

the leading order term of these two formulae match up, but Bogomolny et al. show 
that by looking at a scaled version of the statistic and then choosing an "effective" 
matrix size, related to A^ by multiplication by a constant of arithmetic origin, they 
can match the first lower-order term in the formulae. By scaling a further variable 
they match yet another lower-order term, although this step requires that the con- 
sidered Riemann zeros be high lying, a point we will return to later. The authors 
illustrate with comprehensive numerical results that with each successive refine- 
ment the fit of the random-matrix model to the Riemann zero data significantly 
improves. What is remarkable about this is that they don't just see good numeri- 
cal agreement in the two-point correlation statistic, where good agreement is to be 
expected as terms were matched by design, but also in the nearest-neighbour spac- 
ing distribution. The nearest neighbour spacing statistic is the probability density 
for distances between consecutive zeros, or equivalently, a normalized histogram 
of gaps between consecutive zeros. The argument is that the corrections to the 
asymptotics in the form of A^eff and the further scaling factor is also valid for all 
correlation functions, and therefore also for the nearest neighbour spacing. The 
strength of this work is that using this heuristic approach the information gained 
from a simpler statistic (two-point correlation) yields information for a more com- 
plicated one (nearest-neighbour spacing) — the latter one being determined by all 
correlation functions together. 

In the following we adopt the above method and apply it to the case of J-'^(X). 
In this case we have a conjectural formula for the one- level density for the scaled 
zeros of J-'^(X), as derived in [371|38]. For our purposes this is best given as the 



Random matrix theory does not see the fine arithmetic properties of the family, which surface 
in the lower-order terms. Thus, while the main terms of various families of elliptic curves are 
the same, the lower-order terms show differences due, for instance, to complex multiplication or 
torsion points; see [CTj . 
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following expansion for large X (equation (3.18) in [3^): 



^ ^ 5:h^) 



0<d<X 7d 

Xd(-M)wB=+l 

... sin(27rr) 1 + cos(27rr) vrr sin(27rr) \ , ^/ 1 



with 



(4.2) 



logl^^TT^J' (4-3) 



MX\ 

2^ 



where jd is the imaginary part of a generic zero of LE{s,Xd), the sum is over 
fundamental discriminants d, and X* is the number of fundamental discriminants 
satisfying the conditions on the sum. The coefficients ri and r2 in 04. 2 p are arith- 
metic constants involving the Dirichlet coefficients of Le{s); see [3Z] for details. 

As expected the result (14. 2 p has the form of the asymptotic random-matrix 
result 

7^ , N /'sin27rs\ , , ,. 

for the even orthogonal group, plus correction terms. We compare this expansion 
with the one we obtain by expanding the scaled one-level density of S0(2A^) for 
finite A^. The unsealed one-level density of S0(2A^) is (see, for example, j9j) 

111 It: sm s 

and so scaling by the mean density and expanding in powers of 1/A^ gives 

^ Z3 ('^y\ 1 I sin(27r?/) 1 + cos(27r?/) vry sin(27ry) . ^( ^\ , . „^ 
N'^A-N)=^^^2^y 2iV 6iV^ + ^ ^iV^ J ' ^^'^^ 

By choosing an effective matrix size 

A^eff = ^ (4.7) 

we match the next-to- leading term in (14. 2 p and (14.61) . Arguing as Bogomolny et 
al. we conjecture that the improvement made by using matrices of size Agff also 
holds for all n-point correlation or density functions. In particular we apply this 
result to the distribution of the lowest zero and see significantly better agreement 
in the bulk and tail of the distribution when we use Agff as opposed to the use of 
the standard matrix size A'std = ^■ 

We illustrate the effect of Aefi in figure HI We choose the elliptic curve Ex\ 
((ci, C2, C3, C4, Cg) = (0, —1, 1, 0, 0) in Weierstrass form) and find numerically that, 
for this curve, 

ri ^ 2.8600. (4.8) 
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We used Rubinstein's Icalc [56j to compute the lowest zero of the quadratic twists 
Leu ("S, Xd) with even functional equation and with fundamental discriminants d, 
Q < d < X = 400,000. The standard matrix size corresponding to X = 400,000 is 



iVstd = log 
whereas the effective one is 



IIX 



27r 



12.26, 



K 



N. 



std 



eflf 



2ri 



2.14. 



(4.9) 



(4.10) 



The distribution of the lowest zero (a normalized histogram of the heights of 



1.5 



0.5 




0.5 



1.5 



Figure 4. Distribution of the lowest zero for LEi^{s,Xd) with < 
d < 400,000 (bar chart), distribution of the lowest eigenvalue of 
S0(2X) with Xeff = 2.14 (solid), standard Xstd = 12.26 (dots). 

the lowest zero of each L-function in the family) is then depicted (for rank-zero 
curves) as a bar chart in figure H] whereas the distribution of the lowest eigenvalue 
of S0(2X) with 'effective' matrix size Xgff is the solid curve and the one with 
'standard' matrix size N^td is the dotted curve. The distribution of the lowest 
eigenvalue is related to the solution of a non-linear ordinary differential equation 
of Painleve VI (see [23]). The curves for N^s and N^td in figure [Dare computed 
for these non-integral values of N by establishing a numerical differential equation 
solver for Painleve VI, which is developed in [20j . 

We observe from figure H] that the distribution of the lowest eigenvalue of 
S0(2A^cff) niimics the distribution of the lowest zero of rank-zero curves from 
J-'^(X) much better both in the bulk and in the tail of the distribution than 
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S0(2A^std). However, near the origins we still find a large discrepancy. We address 
this discrepancy in the next section. 
Substituting N^s into fl4.6p we obtain: 

'^ ID f^y\ f, , sm{2TTy)\ l + cos{2iTy) Arf ny sm{2ny) ,^/l\ 

(4.11) 
With the substitution A^ = N^s the scaled one-level density of S0(2A^) now agrees 
with our conjectural answer f l4.2p for that of J-'^(X) in the leading and next-to- 
leading order term. The method of Bogomolny et al. indicates how to proceed 
to find agreement down to the next order term of order L~^. However, with our 
available data for twists up to X = 400,000 we are not working with sufficiently 
large X to proceed further with their method. For completeness, the details are 
worked out in 



5. Cut-off value for the excised ensemble 

5.1. Calculating the cut-off value. We now develop an argument to determine 
the cut-off value c in f l3.4p . In jTS] and [S], Conrey, Keating, Rubinstein and 
Snaith developed a method using random matrix theory to conjecture the asymp- 
totic order of the number of L-functions in a family such as J-'^(X) that vanish at 
the central point; alternatively, this is equivalent, on the Birch and Swinnerton- 
Dyer conjecture, to the asymptotic order of the number of curves of rank 2 or 
higher in the associated family of elliptic curves. Numerical tests support the 
prediction for the order, but the associated proportionality constant could not 
be determined. We will use numerical findings in those papers and the method 
introduced there to arrive at the cut-off value c. 

We begin by reviewing the method of [13], modifying it to meet our current 
purpose. Adjusting slightly the notation of [13] to ours, let 

Me{X,s) = ^^ E LE{l/2,Xd)\ (5.1) 

0<d<A' 

i'i5(».xd)e-^i!;(x) 

where X* = ^{0 < d < X \ LE{s,Xd) ^ -^Ei-^)} ^^ ^^^ number of terms in the 
sum above. Following the philosophy set out in ^2] and [13] we expect that, for 
large X and A^ ~ logX, 

Ms(X,s)~a,(£;)Mo(iV,s), (5.2) 

where 



Mo{N,s)= AA{l,NydA (5.3) 

'sO{2N) 
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and as{E) is an arithmetical expression depending on the Dirichlet coefficients 
A(p) of the curve E: 



as{E) 




C 



P\pl/2 



c 



-1 



P\pl/2 



(5.4) 



.p]M- 

This holds for prime conductor M, where uj{E) is the sign of the functional equa- 
tion of Le{s). For our purposes, the ± in the last line above must actually be +, 
corresponding to twists by positive fundamental discriminants < d < X (the — 
sign corresponds to twists by negative fundamental discriminants —X < d < 0). 
We define 



Cp{z) 



J2hp 



"^2" 



\{p)z + i)M{p)z^ 



(5.5) 



with ipuip) given at 02.51) . 

In 05. 2p . Mo{N,s) denotes the moment generating function of the values |A^(1, A^)| 
as A varies in the random matrix ensemble S0(2iV) (the s^^ moment is the ex- 
pected value of |A^(l,iV)|'^). For Re(s) > —1/2, Mo{N,s) can be explicitly eval- 
uated HSl as 



Mo{N,s) 



N 



AA{l,NydA = 2 



2Ns 



SO{2N) 



yr m + J-l)T{s+J-l/2) 

lir(j-i/2)r(. + j + iv-i)- ^ • ^ 



Thus, from 05.61) . we have information about the value distribution of the charac- 
teristic polynomials. More precisely, we have, for c > 0, that 

"C+ioo 



PoiN.x) 



27iix 



Mo{N,s)x~'ds; 



(5.7) 



here Po{N, x) denotes the probability density for values of the characteristic poly- 
nomials Aa{1,N) with A e S0{2N). 

For small x > 0, the regime of our interest, the major contribution in the integral 
on the right side of 05. 7p comes from the simple pole at s = —1/2 in (15. 6p . thus 

Po{N,x)^x-^/^h{N) (5.8) 

where 

h{N) = Res Mo{N, s) = 2-^r(iv)-i f] r(iV + j - l)r(j) 
«=-i/2 ^^ ' ^ ^ ' y r(j - i/2)r(j + A^ - 3/2) 

We also make use, for large A^, of the asymptotic 

h{N) ~ 2-^/«G'(l/2)7r-i/4Ar3/8, 
where G is the Barnes G-function J2j. 



(5.9) 



(5.10) 
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Notice that Po{N,x) dx is the probabihty that a characteristic polynomial 
Aa(1, N) oi a E S0{2N) takes a value between x and x + dx. Hence 

Prob(0< Aa(1,A^) <p) = / Po{N,x)dx. (5.11) 

Jo 

With fIS.Sp we have for small x that 

rp 
Prob(0 < AA(l,iV) <p) ~ / x-^/^h{N)dx = 2p^/^h{N). (5.12) 



We now define another probability density PE^d, x). Due to 05. 2p . we expect that 
this is in some sense a smooth approximation to the probability density for elliptic 
curve L- values from J-'^(X) which have fundamental discriminants around d. We 
define, in analogy with (15. 7p . 

PE{d,x) := -— / as{E)Mo{\ogd,s)x-'ds ~ a_i/2{E)Po{\ogd,x). (5.13) 

ZTXIX J c—ioo 

The final approximation above holds for small x; it is obtained by shifting the line 
of integration left past the pole at s = — 1/2 and picking up the respective residue 
of the integrand. 

As described in section [3l the formula of Waldspurger et al. (12.70 implies a 
discretization of central L- values for elements from J-'^(X) given by 

L£;(l/2,Xd) = whenever L£;(l/2, Xd) < -^. (5.14) 

By calculating the probability that a random variable, Y^, with probability density 
PE{d, x), takes a value less than KE/\fd and summing over d up to X, the authors 
in [13] and [H] predicted the correct order of magnitude for the number of L- 
functions that vanish at the central point (see Conjecture 5.1 of [H]). However, 
the constant factor could not be predicted correctly to agree with the numerical 
data. Since their cut-off value KE/yd did not give the correct number of vanishing 
L-functions, we will presently work backwards and use the numerically calculated 
number of L-functions that are zero at the central point (numerical data from [13]) 
to deduce an "effective" cut-off -jB- that does give the correct answer, for some 
5 > 0. Thus we repeat the calculation of [13] but with the modified cut-off value. 
We have 



< Id < —j^ I ~ / a^xi2{E)x-^l'^h{^ogd)dx 

/r \ 1/2 

= 2a_i/2(i5) (^^j h{\ogd). 



1/2 (5-15) 
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Now we follow fT3] and conjecture that 

MLe{s, Xd) e J-+(X), d prime : Le(1/2, Xd) = 0} = $^* Prob (o<Y,< ^"j 



d<X 
d prime 



-^5^2a_,/,(i?)(^^j MlogX), (5.16) 



where the starred sum means that d is restricted to prime fundamental dis- 
criminants for which Xd{—M)uj{E) = +1 (of which, asymptotically, there are 
X/(41ogX) of size at most X). Using (IS.lOp . we get 

#{Ls(s, Xd) e J-+(X), d prime : L£;(l/2, Xd) = 0} 

~ jj^ ■ 2 a_i/2(i^)v^ 2-^/«G(l/2)7r-V^(logX)^/«<5i/2 ■ fX^/^ (5.17) 

We wish to obtain a numerical value for 6. We would like to thank Michael 
Rubinstein for sharing his data from [H] on the number of L-functions that vanish 
at the central point. For a large number of elliptic curve families he computes the 
left side of (I5.17P and divides it by 

^a^,/2{E)^^X'/\\ogXr'/\ (5.18) 

The results are plotted (see [11]) as a function of X; the curves flatten out and 
seem to approach a constant value. One such constant is: 

0.2834620 for E = UAr (5.19) 

where the nomenclature of the elliptic curve being twisted to form the family refers 
to Table 3 of [H]. We will refer to this family as twists of En. 
Thus for twists of Eu we have 

|2"^/8^(l/2)7r-^/^5^/' ^ 0.2834620. (5.20) 

With G(l/2) evaluated as approximately 0.603244, we have 

5^0.185116. (5.21) 

For this family we have 

ke = 6.346046521 and a_i/2(^) = 0.732728078. (5.22) 

Thus, with 6 given in (I5.2ip and ke in (I5.22p . we take the 'effective' cut-off to be 

S ■ ke = 1.17A75. (5.23) 

That is, when modeling the distribution of L- values with PE{d,x), as described 
in [13], integrating up to the value of 6KE/Vd gives the correct number of L- 
functions taking the value zero at the central point. 
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5.2. Numerical evidence for the cut-off value: Standard A^std- The proba- 
bility densities Pe and Pq are related, for small values of x (see 05.81] and (15.131) ). 
by 

PEid,x) ~ Po (iV,td,a:2/2(^)x) (5.24) 

when iVstd ~ logd. Thus a cut-off of 6KE/vd applied to PE{d,x) scales to 

Ctd X exp(-Ar,td/2) := aZl/^{E) 6 ke ^ exp(-A/'std/2) (5.25) 

for the distribution of values of characteristic polynomials of matrices of size A'^std- 
Substituting the numerical values we obtain 

Cstd ~ 2.188. (5.26) 

We now present data using the standard matrix size A'^std to model the zero 
statistics. Although we expect that the excised orthogonal ensemble models L- 
functions with discriminant around the value X, for numerical tests we take all d < 
X and set A'^std ~ ^ogX so that we have a substantial data set. 

As mentioned in section HJ and illustrated in figure HI eigenvalues of S0{2N) 
matrices with A^ ~ logX do not give particularly good agreement with the sta- 
tistics of zeros of L- functions from J-'^(X). The procedure of excising matrices 
with small values of |Aa(1,A^)| from the ensemble does not substantially change 
the bulk of the distribution. However, for En, we illustrate in figure |5] how we can 
scale the mean of the distribution of the first eigenvalue to obtain better agree- 
ment. The distribution of the first zero of even quadratic twists of Le^^^s) by 
prime discriminants between and 400,000 has a mean of 0.4081. The distribu- 
tion of the first eigenvalue of 3 x 10^ matrices from S0(2A'std) conditioned to have 
AA(l,Ar,td) > 2.188 X exp(-Ar,td/2), with A^.^d = 12 (^ log(v^ 400,000/27r)) 
was found to be 0.365; deviation from the distribution of the first critical zero of 
the quadratic twists is quite visible, as illustrated in the left graphic of figure O 
Rescaling the random-matrix mean to match that of the L-function zeros gives 
the graphic on the right side of figure |5l showing much better agreement. The 
cumulative plot of the scaled distribution is shown in figure [61 

Working with the mean-scaled RMT results for A^ = A^std, we test various values 
of c in (13. 4p . We measure the agreement with the zero distribution by averaging the 
absolute value of the difference between the cumulative distribution of the zeros 
and the cumulative distribution of the eigenvalues at a set of evenly spaced points 
(the blue crosses on figure E]). The plot of the difference between the cumulative 
distributions versus the cut-off parameter is shown in figure [TJ where c varies along 
the horizontal axis. We see a minimum at Cstd = 2.188; this is the value predicted 
in section [STTj and is marked with a dotted vertical line. 

Initially we thought that an equally plausible way to calculate a cut-off value 
would be to equate not the probability densities, as in (I5.24p . but rather the prob- 
ability of finding an L-value of zero with the probability of discarding a random 
matrix on the basis of the condition (13. 4p . Equating probabilities would lead to a 
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Figure 5. Probability density of the first eigenvalue from 3 x 10^ 
numerically generated matrices A G S'0(2A^std) with |A^(1, A^std)| > 
2.188 X exp(— iVstd/2) and iVgtd = 12 (red dots) compared with the 
first zero of even quadratic twists Lej^^ {s, Xd) with prime fundamen- 
tal discriminants < d < 400,000 (blue crosses). In the left picture 
the random matrix data is not scaled, in the picture on the right 
the mean of the distribution is scaled to match that of the zero data. 




Figure 6. Cumulative probability density of the first eigenvalue 
from 3 X 10^ numerically generated matrices A G S'O (2 A^gtd) with 
|Aa(1, iVstd)| > 2.188 X exp(-iVstd/2) and A^'^td = 12 (red dots) com- 
pared with the first zero of even quadratic twists Le^-^ (s, Xd) with 
prime fundamental discriminants < d < 400,000 (blue crosses). 
The random matrix data is scaled so that the means of the two 
distributions agree. 



calculation such as 

"SKs/Vd 



PE{d,x)dx 



Ske/v^ 



a^i/2{E)x ' h{\ogd)dx 

a^_i^^{E)5KE/Vd 

y~''^h{N,,^) dy 



(5.27) 



o?i-^,2{E)S>^E/Vd 



Po{N^td,y)dy, 
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Figure 7. A measure of the difference between the mean-scaled 
cumulative distribution of the first eigenvalue of the excised random 
matrix model f l5.26p with A^std = 12 for various values of c and the 
cumulative distribution of the first zero of even quadratic twists 
LEii{s,Xd) with prime fundamental discriminants < d < 400,000. 
The values Cgtd = 2.188 is marked with a vertical line. 



implying a cut-off of 

a\/^{E)SKE eM-Nstd/2) ^ 0.6307 exp(-iV3td/2). (5.28) 

From figure[7]we see that c = 0.6307 certainly does not minimize the error between 
the excised random matrix model and the zero statistics. 

It is noteworthy that matching probability densities, rather than probabilities, 
appears to be the correct model. Taking the cut-off value Cgta that best models 
the zero data (as calculated at the beginning of this section by matching densities) 
means that the proportion of the matrices that are being excluded by the condition 
(13. 4p is not the same as the proportion of L-functions from our family of quadratic 
twists that are excluded from the zero statistics because they have rank higher than 
zero. The reason for this, and whether the excluded matrices might potentially 
model the L-functions associated with the higher-rank curves, are topics for future 
investigation. 

5.3. Numerical evidence for the cut-off value: A'efr. The excised model can 
be applied using orthogonal matrices of any even size 2N. We will presently use a 
method identical to that in the previous section; however, now we find the cut-off 
for matrices of size N = N^ff (which was calculated in section H]). We recall the 
shape of Po{N,x) for small x and large N, equations (15. 8p and (I5.10p . and see 



PE{d,x) ~ a^i/2X-^/^h{\ogd) ~ a^i/2X-^^^h{2riN,s) 



~ a^y2X-'^\2nf'h{N,s) 



(5.29) 



P< 



o 



N,s,aZl/2iE){2n)-'/^x 



Thus a cut-off of 5KE/\/d applied to PE{d,x) scales to 

Ceff X exp(-Arstd/2) := aZl/2{E){2ri)-^/^ 6 ke x exp{-riN^s) 



(5.30) 
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for the distribution of values of characteristic polynomials of matrices of size A'eff . 
Plugging in the numerical values we obtain, with ri = 2.8600 for Eu (given at 

(W) 

Ceff ^ 0.5916. (5.31) 

For this matrix size we use the cut-off value c = Ceg in (13. 4p . 

We now compare data for zeros of elliptic curve L-functions to matrices of 
size A^eff- 

We note here that, although the most accurate description of our expectation is 
that the excised orthogonal ensemble models L-functions with discriminant around 
the value d when N^g ~ logd/{2ri), for numerical tests we take all d in an interval 
such as < d < X and set Nes ~ logX/(2ri). This is the most effective way to 
acquire enough data to have good resolution in the plots. 

For Ell we have computed zeros of even twists L^;^^ (s, Xd) for prime d between 
and 400,000 which do not vanish at the central point. We compare the distribution 
of the first zero above the critical point of these L-functions with the numerically 
generated distribution of first eigenvalues of 3 x 10® matrices from the excised 
ensemble S0(4) (in section H] we computed A^eff = 2.14) for various values of the 
cut-off c. Figure [H] shows the probability density and the cumulative probability 
density of first zeros and eigenvalues for Ceg = 0.5916. In contrast to modeling with 
matrices of size iVstd, here we do not scale the mean of the distribution; however, 
for those interested, we note that the probability density of the zeros has mean 
0.4081 and the probability density of the eigenvalues has mean 0.4234. 




Figure 8. Probability density (left) and cumulative probability 
density (right) of the first eigenvalue from 3 x 10® numerically 
generated matrices A e S0{2Nes) with \AA{l,Nes)\ > 0.5916 x 
exp(— riA^cff) and A^eff = 2 (red dots) compared with the first zero of 
rank zero even quadratic twists Le^ {s, Xd) with prime fundamental 
discriminants < d < 400,000 (blue crosses). 

In figure |9] we plot a measure of the difference between the first zero distribution 
and the distribution of the first eigenvalue for various values of the cut-off c. Here 
the error is calculated by summing the absolute value of the difference between 
the cumulative distribution of the zeros and the cumulative distribution of the 
eigenvalues at a set of evenly spaced points (the positions of the blue crosses in 
figure E]) and dividing by the number of points. We see a minimum at Ceff = 0.5916; 
this is the value predicted in section [5?T] and is marked with a dotted vertical line. 
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Figure 9. A measure of the difference between the cumulative dis- 
tribution of the first eigenvalue of the excised random matrix model 
(I3.4p . as c varies along the horizontal axis, and the cumulative distri- 
bution of the first zero of rank zero even quadratic twists Le^^ (s, Xd) 
of with prime fundamental discriminants < c? < 400,000. The 
value Ceff = 0.5916 is marked with a vertical line. 



As in the previous section, another possible model is to equate probabilities 
rather than probability densities. For A'^eg this leads to 



(2ri) 



3/4^2 



-1/2 



{E)5ke ^ 2.3328. 



(5.32) 



Again, from figure Owe see that this does not come close to minimizing the error. 

5.4. Excised model and the one- level density. In section H] we provided evi- 
dence that the lower-order terms of the one- level density of J-'^(X) determine an 
effective matrix size N^s such that the distribution of first eigenvalues of S0(2A'efj) 
models the bulk and tail of the distribution of first zeros of J-'^(X). In sections IHTTI - 
15.31 we discussed how to obtain an appropriate cut-off value for the characteristic 
polynomial A^ at 1 with A G S0(2A^) so that the distribution of first eigenvalues 
of this subset of S0(2A^) models the region at and near the origin of the distribu- 
tion of first zeros of J-'^(X). Here the matrix size is either the effective one, i.e., 
N = Nes or the standard one, i.e., A^ = Ngtd- In the first case no further scaling 
is undertaken whereas in the latter one we match mean densities of S0(2Astd) to 
J^e{^) to achieve qualitative and quantitative agreement throughout the origin, 
bulk and tail between the cumulative distribution of first eigenvalues of S0(2A^) 
and cumulative distribution of first zeros of J-'g(X). In this section we explore 
to what extent the values for the matrix size and the cut-off value on the char- 
acteristic polynomial used for the distribution of first eigenvalues (to model the 
distribution of first zeros) can be applied to the one-level density of eigenvalues 
(to model the one- level density of zeros). We compare zero data of the fam- 
ily J-'^(X) of quadratic twists LEj^j^{s,Xd) with prime fundamental discriminants 
< d < X = 400,000 to eigenvalue statistics of random matrices from S0(2A^) 
for two values of N. Figure [10] depicts as a red solid line the one-level density of 
matrices A from S0(2Astd) with characteristic polynomials constrained to obey 
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1 -level density: excized RMT, 24x24, cut-off=0. 005424 

1 -level density: excized RMT, 4x4, cut-off=0.001466 

1 -level density: NT, prime conductor 



Figure 10. Red solid line: One-level density of 9.12 x 10^ numeri- 
cally generated matrices from S0(2A^std) and cut-off |Aa(1, A^std)| > 
2.188 X exp(-A/'std/2) = 0.005424, scaled so that the mean of the 
first eigenvalue matches that of the first zero. Black dotted line: 
One-level density of 9.8 x 10^ numerically generated matrices from 
S0(2A/'eff) and cut-off |Aa(1, iVeff)| > 0.5916 x exp(-iVstd/2) = 
0.001466. Blue crosses: One-level density for zeros of even quadratic 
twists LEii{s.,Xd) with prime fundamental discriminant between 
and 400,000. 



|AA(l,A/'std)| > 2.188 X exp(-iVstd/2) = 0.005424 with N^td = 12. We obtained 
this by generating 9.12 x 10^ matrices from SO(24). This one-level density on the 
random matrix theory (RMT) side has been scaled, in the horizontal direction, 
with the scaling factor 0.4081/0.365 = 1.118 which we obtained numerically when 
comparing with the distribution of the first zeros in section 15.21 The one-level 
density of matrices A from S0(2iVeff) with A^cff = 2 where the characteristic poly- 
nomials are constrained to obey |Aa(1, iVes)! > 0.5916 x exp(— A^std/2) = 0.001466 
is depicted as a dotted line. We obtained this by generating 9.8 x 10^ matrices 
from SO (4). No scaling of means has been undertaken in this case. The one- level 
density for zeros of rank- zero, even quadratic twists of En with prime discrimi- 
nant between and 400 000 is depicted as blue crosses. In the A^std case we find 
correlation with the zero data over a wide range whereas in the Nes case we only 
find some agreement up to the first unit mean spacing. Thereafter we see a clear 
discrepancy. 

In figure [11] we consider the associated cumulative distributions from figure [TOl 
Here we observe that the distributions of A^std (red line) and the zero data (blue 
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cumulative distr: RMT, 24x24 

cumulative distr: RMT, 4x4 

cumulative distr: NT, prime conductor 




Figure 11. Cumulative distributions of one-level densities from 
figure [TOl Red solid line: numerically generated matrices from 
S0(2iVstd), with scaling; dotted line: numerically generated matri- 
ces from S0(2A^eflf); blue crosses: zeros of even quadratic twists 
Leii {s, Xd) with prime fundamental discriminants between and 
400,000. 



crosses) track each other over a wide range, whereas the distributions of A'^efr 
(dotted line) and the zero data (blue crosses) track each other only up to the first 
unit spacing. Thereafter the behaviour of the zeros is not captured. 
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1.6 



cumulative distr: RMT, 24x24 

cumulative distr: RMT, 4x4 

cumulative distr: NT, prime conductor 




Figure 12. Magnification of figure [TTl i.e., cumulative distribu- 
tions of one-level densities from figure [TOl Red solid line: numeri- 
cally generated matrices from S0(2iVstd), with scaling; dotted line: 
numerically generated matrices from S0(2A'"eflf); blue crosses: zeros 
of even quadratic twists L^;^^ (s, Xd) with prime conductor between 
and 400,000. 



Figure [T2] is a magnification of figure [TTl Here we observe that also the slope of 
the zero data near the origin is nicely captured by the excised RMT model with 
standard matrix size iVstd- This feature is not so clear in the N^s case. 

In summary, to answer the question to what extent the values for the matrix 
size and the cut-off value of the characteristic polynomial used to model the dis- 
tribution of first zeros by the distribution of first eigenvalues can be applied to 
the one-level density, we find that we achieve better agreement when choosing 
the excised ensemble with standard matrix size (with the data scaled so that the 
mean value of the first eigenvalue matches that of the first zero) than choosing 
the excised ensemble with effective matrix size. In the latter case we only obtain 
agreement up to the first unit mean spacing and in the former we get agreement 
over a wide range. It should be noted here that we only require our model to 
give useful predictions over a distance of one mean spacing, because further from 
the origin we know from [37] that the one level density is strongly dominated by 
arithmetic contributions that are accurately modelled by methods directly incor- 
porating number theoretical information. 
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6. One-level density for the excised random matrix model 

In the previous sections we provided evidence that we can use eigenvalue statis- 
tics of random matrices from the excised orthogonal ensemble to model the zero 
statistics for the family J^^(X) of even quadratic twists of a fixed elliptic curve 
L-function. Our goal in this section is to compute the one-level density Ri^ for 
Tx, the ensemble defined in the introduction. The main results are contained in 
the three theorems set out in the introduction. We prove those theorems in this 
section. 

6.1. Proof of Theorem II. IL 

Proof. Consider 



R^''{9,):=C;,-N r ■■■ r H{\og\AA{l,N)\-X)x 
Jo Jo 



X YlicosOj - cosekyde2- ■ ■ dON, (6.1) 
j<k 

which is the one-level density for the set T^- Here H{x) denotes the Heaviside 
function 

^^, , I 1 for X > , , 

H{x) = 1 (6.2) 

I tor X < 0, 

and Cx is a normalization constant (which we discuss later). Next, we replace 
in 06.11) the Heaviside function with its integral representation 

ipexpM^^^ (6.3) 

where c > 0, and observe that 

TV 



AA{l,N) = AA{exp{te),N) 



0=0 



n(l-exp(^0,))(l-exp(-z0,)) 



i=i 



(6.4) 



N 

Thus, we have 

if (log |A^(1, N)\-X) = ^ r^^ 2^,exp(-r^ ^ _ ^os O^Ydr. (6.5) 

Substituting (16. 5 p into (16. ip gives 

27rz Jc-ioo r Jo Jo j-J^ 

X JJ(cos 9j - cos 9k)^de2 ■■■dON dr. (6.6) 

3<k 
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Now, observe that in (16. 6p we have 
^ " Yl{l-cos9jYY[{cos9j-cos9kyde2---deM 

/•TT /"TT ^^ 



where w*^"''^-'(cos^) = (1 — cos0)""'"^/^(l + cos^)^"*"^/^ is the weight function for the 
Jacobi ensemble of random matrices ^24j . We now observe that 

"TT ...fr-l/2.-l/2')/„„„ /I \ TT/„„„ /I „„./l\2j/i j/^ 

(6.8) 
is the one- level density for the Jacobi ensemble J^. This completes the proof. D 






6.2. Computation of the normalization constant Cx- Recall that when in- 
tegrating over S0(2A^) with respect to the normalized Haar measure using Weyl's 
integration formula [61j the normalization constant Cso{2N) is determined by 



1=1 dA = Cso{2N) I Y\{cos Oj - cos OkYdOi ■ ■ ■ dON- (6.9) 

JSO(27V) J[0,-k]N .^j^ 

Selberg's integral formula states that for integral A^ and complex r, s with Re(r), 
Re(s) > -1/2 

/ d(f)i--- (i0Ar]^(l -cos 0/)'(l + cose/);)'' J]^ (cOS 0j - COS 0fc)^ 

•^^ -^O 1=1 l<j<k<N 

_ N(N+r+s-l) .. n ^^^ + ^^^^' + ^/^ + ^^^^'' + ^/^ + ^^ r6 10) 

Using (I6.10p the normalization constant has the explicit form 

Likewise, the normalization constant Cx is determined by 

l = Cx f H{\ogKA{l,N)-X)W{cosej-cosekfdei---deN. (6.12) 

Using (16. 5 p yields 

Cso^^Csoim f de^---de^dax 

Cx 2Tli Jc-ioo J[0,7r]N 

^ ^Na ^M-O^^) -Q(^ _ cosgj)" J](cOsg, - COsgfc)l (6.13) 
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Substituting (16. lip into ( I6.13P and applying Selberg's integral formula f l6.10p we 
obtain 

Csoi2N) ^ J_ r^'^ eM-aXl,^^ ^' r(iV + j)r(a + 1/2 + j) 

C;. 2ni A_,^ a 11 r(a + AT + j)r(l/2 + j) " ^ " ^ 

The evaluation of fl6.14p boils down to closing the contour to the left and to 
computing the residues associated to the poles at a = and the negative half 
integers coming from the T{a + 1/2 + j)-term. The residue at the simple pole 
a = is 1 and the residue at the simple pole at a = —1/2 is 

2 cxnf.r /212~^ ^(^^ n r(iV + j)r(j) 

-2exp(A'/2)2 r(iV - l/2)r(l/2) H V{N + j - 1/2)1(1/2 + j) " ^^"^^^ 

We denote the contribution to (I6.14p from residues of higher order poles at a^ = 
=^ with A; = 1,2,3,... by 

J]afcexp((A; + l/2)x). (6.16) 

fe>i 

Thus, we obtain 

Cso(2^)_, exp(A'/2) r(iv) ^' r(iv + j)r(j) 



n 



Cx 2^-1 T{N - l/2)r(l/2) J-^-L T{N + j - 1/2)1(1/2 + j) 

+ 5^afcexp((A; + l/2)x). 
fe>i 

(6.17) 

Notice that when X — )• — oo we have Cx — ?■ C'so(2Ar)5 as expected. The regime 
of interest for us is when A* < 0. From (I6.16P we see that the contributions of 
residues of higher poles decrease exponentially for A* < 0. The computation of 
residues of higher order poles is easily done with a computer algebra system (we 
used Mathematica) . We find that for reasonable X we have good convergence 
using about 10 poles. When using only the K rightmost poles the error term is 
0(exp(— Ci^A")) with ax+i < ck < dx- We refer to section |6^ for the relevant 
details. There we will also discuss the convergence of a series similar to (I6.16p . 
We are then left to show that the contribution from closing the contour is 0. This 
is indeed the case for exp(A') < 2^^ = max[A^(l, A^)] but we skip the calculation 
here and instead also refer to section WM as we do a similar computation for i?^-* 
there. 

6.3. Proof of Theorem 11.21 

Proof. We rewrite (16. 6 p so we can apply standard methods in the theory of or- 
thogonal polynomials. First we consider the general Jacobi ensemble with weight 
function (1 — x)°'(l + x)^ and then specialize to our setting with a = r — 1/2 and 
/3 = —1/2. Notice that the weight function here differs slightly from the one in 
section [TJ 
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Following Szego [59] we define the Jacobi polynomials Pn (x), a, /3 > — 1, 
which satisfy 

Pt'^\x)PJ^'^\x){l - xr{\ + xfdx = 5nmhl:'^\ (6.18) 



where 



/zl"'« 



2«+/3+i Y{n + a + l)r{n + P + l) 



2n + a + (3 + lT{n + l)T{n + a + /3 + l)' 
Changing to angular variables, f l6.18p takes the equivalent form 



(6.19) 



P^'''^\cOSe)PJ^''^\cOSe){l-COSe)''+^/\l+COSef+^/^de = 6nmhl^'^\ (6.20) 



Using the Vandermonde determinant and further matrix row and column oper- 
ations we have 



1 |(cos6'fc — cos^j 



j<k 



1 cos 01 cos^ 6i 
1 cos 62 cos^ 62 

1 cos 9n cos^ 9]\r 



cos^-i ^9 



cos^-i 9 



N 



N-1 f7^ia,l3)\l/2 

^ TT^-^ — 

i=o S 



X 



>{",/9)/ 



,{",/3) 



^Pr''^^(cos^i 



,(",/3) 



,(<^,/9) 



PArTAcos^i) 



,{",/3)/ 



=^Po^"''^^ (cos ^2) ^^P{"''^^ (cos 92) 



,(«,/3)/ 



bT^O^"'^^(-S^.) ^rfe^Pr^(cOS^ 



(«,/?)/ 



u(c,/3) 



('("./9) ; 



^ 



c,/3) 



7Af ; 



1 u('^^P)t^^r.a 



where i) is the leading coefficient of the polynomial P, 



(",/?). 



] 



^(a,/3) ^ 2-i Z^^-^' + « + /^ 



(6.21) 



(6.22) 



Thus we have the following determinantal expression 

I [(cOS^fc — COS^j)^ 



j<k 



N-1 T(a,/3) 



n 



h\ 



N 



Setting 



{"./9)\2 
J > 



det 5:(e?)-^^i-f (cos^,)Pi"_f (COS0, 



vn=l 



j,fc=l,...,Af 



Af-1 



^ .^ 2-^(^+''+«-i) TT 



V{s + r + N + i] 



j=0 



T{2 + j)T{s + 1/2 + j)T{r + 1/2 + j] 



(6.23) 



(6.24) 
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the normalized measure is 

Cs,r n(l - cos%)'"(l + cosOjY YlicosOk - cosOj)^ 

j=l j<k 

1 ^ 

= ]^det ((5:(e//^'^-^/^v^e-//^'^-^/^Hcos^,oe-i^/^'^-^/^Hcos^,))x 

n=l 

X (1 + cose,y/\i - cose.y/^i + cosOky^^i - cosOkY'^ 



^,^ei{ft'^^^-'"\0,A 



j,k=l,...,N 



k) I , 

j,k=l,...,N 

(6.25) 



(.(r-l/2,s-l/2). 



where fj^~ ' ' ^ ' '(6j,6k) is imphcitly defined as the expression inside 'det(-)' in 
the middle term above. 

Observe that this determinantal kernel /at satisfies the hypotheses of Gaudin's 
Lemma (see e.g. Theorem 5.2.1 in |16]), namely 

f^^'^'^'''^'^\e,e)de = N (6.26) 

and 

^ fN{x,e)fN{e,y)de = fN{x,y). (6.27) 



'0 
By Gaudin's Lemma we then have 



det(/(^„ eu))j,k=i^...,NMN = {N-{N- 1)) det(/(0,, eu)),,k=i^...,N-i. (6.28) 

Applying (I6.28P N — n times, together with (I6.25p . gives the following formula for 
the n-level density: 



/•TT PIT 1 

Rn{6l, ■ ■ ■ ,0n) = J-rz Tj- / ■ ■ ■ / Try det(/Ar(6'j, Ok))NxNdGn+i ■ ■ ■ (19^ 



1 



(AT - (iV - 1)){N - (iV - 2)) ■ ■ ■ (AT - (n + 1 - 1)) det(/^(^„ 4)), 



{N — n)\ 

= det{fNi0j,9k))nxn- 

(6.29) 
So, using (16.291) in (16. 6p . we arrive at 

Rl^{9) = ^ r^^'-^^^^^ft'^''-'^'\0,0)dr (6.30) 



2vri J^_i^ Co,r r 

with Co,r given in (I6.24p . 
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With the Christoffel-Darboux formula (see equation (4.5.2) in [59]) we have 



I 



(r-l/2,.-l/2)j 



^j,dk) 



2-r-s+l 



T{N + l)T{N + r + s) 



2{N -l) + r + s + l r(A^ + r - l/2)r(A^ + s- 1/2) 



X 



;i + cosejy/^{i - cosOjY/^i + cos 0^)^/2(1 _ cosOkY 



/2 



COS 9j — cos 9k 



P. 



(r-l/2,s-l/2) 



N 



' COS ej)p^_^ 



(r-l/2,s-l/2). 



P, 



COS 9 k) 

(-l/2,.-l/2),^^^,.p(.-l/2,s-l/2).^^. 



. (6.31) 



Thus f{9,9) in (I6.30p reduces, with the Christoffel-Darboux formula (I6.3ip . to 

2^-''' r(A^ + l)r(A^ + r) 



.(r-1/2,-1/2) 
J N 



cos( 



2A^ + r - 1 r(A^ + r - l/2)r(A^ - 1/2) 



X 



d 



dcos9 



P 



(r-l/2,-1/2) 



N 



{cos 9) 



Pt'^''-'^'\^osi 



Ptl^''-'^'\cosi 



d 



acost^ 



)1— r 



where we define 



P{N,r,9) 



cost 



d 



T{N + l)T{N + r) 



2N + r-l T{N + r - l/2)r(A^ - 1/2) 



PiN,r,9), (6.32) 



p(r-l/2,-l/2) 

dcos9 ^ 



(cos( 



P: 



(r-l/2,-1/2), 



N~l 



COS I 



_p(;-i/2,-i/2)^^^^, 



d 



dcos9 



:P, 



(r-l/2,-1/2), 



AT-l 



cos^). (6.33) 



With fl6:32|) in flOOj) we arrive at 

c+joo 



^''■<^'-i 



exp(-rA') ^^2^2JVr-JV.. 



AT-l 

j=0 



X (1-cos^)^ 



r(2 + j)r(i/2 + j)r(r + i/2 + j) 

T{r + N + j) 

2^-^ r(A^ + l)r(iV + r) 

2Ar + r - 1 T{N + r - l/2)r(A^ - 1/2) 



X 



P{N,r,9)dr. (6.34) 



D 



6.4. Proof of Theorem 11.31 

Proof. We consider (16.341) (that is the form of R^'^ given in Theorem ll.2p . First, we 

observe that by [59], pages 63-64, the Jacobi polynomial P^' (x) is a polynomial 
in a, P and x for arbitrary complex values of a and /3. Hence P(iV, r, 9) is a 
polynomial in r. 
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The poles arising from r(A^ + r) at negative integers < —N are cancelled by 
the zeros of the term l/r(r + N + j) for j = 0. The pole at r = 1 — 2N of 
1/[2N + r — 1) is cancelled by the zero of l/r{N + r+j) when j = N — 1. We now 
discuss under which conditions we can close the contour to the left. For this, it 
is helpful to analyze the integrand of fl6.34p when \r\ is large. Our main reference 
for various identities and formulae in the following is [T]. 

Using equation (22.5.42) in [1] we write the Jacobi polynomials as 

pi^-i/2,-1/2) ^^^^ 0)=f^ + y- V2\ F{-N,N + r;r + 1/2; '-f^), (6.35) 

where F is a hypergeometric function. This representation allows us to calculate 
the derivative of the Jacobi polynomial that appears in P{N,r,6) using (15.2.1 
ind]) 

—F(a, b; c; z) = —Fia + 1, 6 + 1; c + 1; 2). (6.36) 

dz c 

Combining the last identity with (I6.35P gives 



d 



,(r-l/2-l/2). 



p;^; ' '" ' ^(cos^) 

■ - 1 /9\ d 

F{-N,N + r;r + 1/2; z] 



dcosO^^ 

1 /N + r- l/2\ d 

~~2\ N )Tz 



1— cos( 



,(N + r- l/2^ (-N)(N + .)^(_^ + 1. ^ + , + 1. , + 3/2^ .-e. 



2V N J r + 1/2 



Substituting (ESSD and dOT]) into ([S33D yields 

P(N rf,)^ ^/iV + .-im/iV + r-3/2 

^ ' ' ' 2(r+l/2)V N )\ N 



2 



(6.37) 



X 



(-A^)(A^ + r)F{-N + 1, A^ - 1 + r; r + 1/2; i^f^) x 

X F{-N + l,A^ + r + l;r + 3/2; i^f^) 
- {-N + 1)(A^ - 1 + r)F{-N, N + r;r + 1/2; i^f^) x 



1— cos6 



X F(-Ar + 2, AT + r;r + 3/2, ^ 

(6.38) 

We want the asymptotics for large |r|, and this is easier when r appears only 
in one argument of the hypergeometric function. For this, we use identity 15.3.4 
from [1]: 

F(a, h; c; z) = {1 - z)-^F{a, c-b;c;j^), (6.39) 
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and so we obtain 



X 



-N){N + r)F{-N +1,-N + 3/2; r + 1/2; 



cosg— 1 ' 

cose+1' 

COS 6—1 N 



X 



X F(-iV + 1, 1/2 -iV;r + 3/2; _,^^, 
- {-N + 1)(A^ - 1 + r)F{-N, 1/2 - N;r + 1/2; ""'^^"^^ 



cos6»+l 



X 



). cos 0—1 ' 



X F(-iV + 2, 3/2 -iV;r + 3/2; _,^^, 

(6.40) 



Recalling equation 15.7.1 from [T] 
F{a, b; c; z] 



E '':^.:!^;^"'y g+ 0(1^1-°-'). 



n! 



(6.4i: 



^^^ T{a)T{b)T{c + n) 

which holds for fixed a,b and z and large \c\ we get 

ir(_Ar + l,-iV + 3/2;r + l/2;£2i|_i)=l + 0(|rri)^ 

and similarly for the other hypergeometric functions. Concentrating on the inte- 
grand of Ri'^{9), in (16.341) . and neglecting factors not depending on r, we find, 
using (16.421) , that the growth-dependence of the integrand on r is 



(6.42) 



22^71- cos ( 



N-l 

n 

i=o 



r(r + l/2+j) exp(-rA') 
T{r + N + j) r 



T{N + r) 



2N + r- ir(A^ + r- 1/2) 



X 



r(A^ + r + l/2)r(A^ + r - 1/2) 



l + 0{\r\-^) . (6.43) 



r(r + l/2)r(r- 1/2) 

Upon simplification, and replacing the products of Gamma functions with Barnes 
double Gamma functions, the growth-dependence on large \r\ is 

exp(r(log(l - cosb") - X)) — +r-. — ._. — -^r-, x 



r{2N + r - 1] 



Gir + 3/2) 
T{N + r + 1/2) 



G{r + 2N) 
l + 0{\r\-^)\. (6.44) 



r(r - 1/2) 

Recall the asymptotic formula for G{z) for large \z\: 

\ogG{z + 1) ~ zW \ogz -\) + \z\og{2Ti) -^^\ogz + C'(-l) + O(z-i). (6.45) 

Applying (I6.45P and Stirling's formula we conclude that the integrand in (I6.34p 
grows like 



^exp(r(2Ariog2 - log 2 + log(l - cos9) - X)) 1 + 0{\r\-^) 



(6.46) 



for large \r\. 

From (16.461) we deduce that if {2N — 1) log 2 + log(l — cos 6) — X>{) then we can 
close the contour in the left half of the complex plane, thus enclosing the poles at 
zero and the negative half integers. If, on the other hand, {2N — 1) log 2 + log(l — 



RANDOM MATRIX MODEL FOR ELLIPTIC CURVES WITH FINITE CONDUCTOR 33 

COS 6) — X < then we must close the contour in the right half-plane, implying 
that the integral f l6.34p is zero. 

Finally, we compute the residues when moving the contour to the left. There is 
a simple pole at r = with residue 

Ai r(iv + j) 2A^-ir(iv-i/2)r(iv-i/2) 

xP{N,0,e). (6.47) 

Note that when r = the expression for P{N, r, 6) reduces to that corresponding 
to S0{2N). 

There is also a simple pole at r = —1/2 with residue 



-2exp(A'/2)2 



N-.2N+3/2 ji-cosey'/' r(iv + i)r(i/2) 



2N-3/2 T{N -1)T{N -1/2) 



N-l 



r(2 + j)r(i/2 + j)r(j) 



xn^^^^(^^^^^i7t^W -1/2.0). (6.48) 

Similarly to f l6.17p . where we determined the normalization constant C^, the 
contribution from the residues of the higher order poles at ak = — ^y^ with 
A; = 1, 2, 3, ... is 

oo 

J2bkexp{{k + 1/2)X) (6.49) 

fc>i 
where the coefficients bk come from the residues. D 

To show convergence we may write the series (16.491) as 

exp{-rX)f{r)dr (6.50) 

where exp(— rA')/(r) is the integrand in (I6.34p and ax+i < ck < c^k- By (I6.46p . 
f{r) -C exp(— rX) for large \r\ and A" < 0. Thus, 

/•Cx+^OO /"OO 

/ exp{-rX)f{r)dr = exp{-CKX)i exp{-itX)f(cK + it)dt 

J cji-ioo J —oo (6.0ij 

= 0(exp(-c^A')). 

Hence, the contributions of residues of higher poles decrease exponentially fast for 
X <{] and using only the K right-most poles gives an error that is 0(exp(— c^frY)). 
The computation of residues of higher order poles can be carried out by computer 
algebra. We find that for reasonable X we have good convergence using about 10 
poles. 

Remark 6.1. The hard gap refiected in Theorem 11.31 may be understood as fol- 
lows. Recalling that Aa(1, N) = 2^ IliLill-cos^n) it follows that logAA(l, A^) = 
(2A^ - 1) log 2 + log(l - cos 6) for the matrix A e S0(2A^) having 2iV - 2 eigenval- 
ues at —1 (i. e., 02, ■ . . ,9j\f = n) and a symmetric pair of eigenvalues at e^*^ (i. e., 
$1 = 6). Plainly, such matrix A maximizes logAA(l,A^) subject to the condition 
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that a pair of eigenvalues lie at e^*^. Therefore, for any fixed X, the condition 
log Aa{1, N)> X implies X < (2N - 1) log 2 + log(l - cos 6): if the latter inequal- 
ity does not hold, there are no matrices A with logA^(l, A^) > X. In particular, 
for fixed N and X there is a lower bound for all eigenphases 9 of any matrix A 
with logAA(l, A^) > A", namely 6 > O^^f := cos-i(l - 2-^'^^'^^^): the excised 
one- level density is identically zero in the excluded spectral interval ( "hard gap" ) 
G < ^inf- Note that the hard gap shrinks to zero exponentially fast as A^ — ?• oo and 
also as A" — !■ — oo. 

6.5. Formula vs. data. For illustration and as a consistency test, in figure [T5] we 
compare our formula for the one-level density with data obtained by generating 
random matrices from the excised ensemble S0(2A^). 



1.2 



0.8 



0.6 



0.4 



0.2 




Figure 13. One-level density of excised S0(2A^),A^ = 2 with 
cut-off |Ayi(l,A^)| > 0.1. The red curve uses our formulae from 
Theorems 11.21 and 11.31 The blue crosses give the empirical one-level 
density of 200,000 numerically generated matrices. 
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